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摘 要 : 基于 


藏 地 区 38 个 气象 站 点 1981 一 2019 年 逐日 气象 资料 ,采用 Penman-Monteith 模型 . 趋 


AHT ~ Morlet 小 波 分 析 、.Mann-Kendall #4 , 4 4 IE X A BC JK (Empirical orthogonal function , EOF) 
和 灰色 模型 探究 潜在 蒸 散 量 (ETo) 的 时 空 演变 规律 以 及 未 来 ET 的 变化 趋势 。 结 果 表 明 : 在 时 间 尺 
度 上 ,西藏 地 区 ET 平均 为 597.12 mm,1981 一 2001 年 呈 显 著 的 减少 趋势 .2002 一 2019 年 呈 显 著 的 
增加 趋势 (P<0.01); 西 藏 全 区 及 五 大 气候 区 年 均 ETo 均 呈现 增加 趋势 ( 除 东 南部 ), 且 均 以 33 a 的 周 
期 震荡 为 第 一 主 周期 。 在 空间 尺度 上 ,年 ETu 主 要 呈现 由 中 部 中 心 向 西南 和 东南 逐步 说 减 的 分 布 
特征 ,高 值 中 心 集中 在 沿江 一 线 , 低 值 中心 位 于 南部 地 区 。 年 ETo 发 生 突变 的 站 点 主要 分 布 在 南部 


边缘 地 区 、 沿 江 一 线 和 东北 部 ,发 生 时 间 集 中 在 201 


Ht 2 80 年代。 构建 的 GM(1,1) 预 测 模型 预测 精 


度 均 值 为 87.85%, 可 用 于 西藏 年 均 ETo 日 期 序列 的 中 长 期 预测 ,预测 结果 显示 , 除 东 南部 年 ETo 的 


预测 值 有 明显 下 降 , 其 他 区 域 均 呈现 上 升 态势 。 


X $ 词 : Penman-Monteith 模型 ; BEAMS; 时 空 分 析 ; 灰色 模型 ; 西藏 


文章 编号 : 


潜在 蒸 散 量 (ET 是 天 气 气候 条 件 决定 的 下 垫 面 
蒸 散 过 程 的 能 力 , 是 土壤 层 水 分 变化 的 重要 指标 ，， 
是 区 域 水 资源 开发 和 评价 依据 ”, 同 时 在 气候 干 湿 
状况 .作物 水 分 僵 亏 和 植被 净 初 级 生产 力 汪 等 方面 
均 得 到 了 广泛 应 用 。 

近年 来 ,诸多 专家 学 者 在 ET 的 模拟 .时空 演变 
规律 .驱动 因素 及 对 气候 变化 的 响应 等 领域 展开 了 
大 量 的 研究 。 秦 鹏 程 等 ”、Wilber 等 拉巴 等 ”、 
钟 巧 等 分别 采 用 Thornthwaite 法 、 人 工 神经 网 络 模 
型 (Artificial neural network , ANN ) .地 表 能 量 平衡 模 
型 (Surface energy balance system, SEBS) , 52 & fs 
(Penman-Monteith, P-M ) ,估算 了 湖北 BUS: Fey Hak 
北 那 曲 、 博 斯 腾 湖 流域 等 地 的 ET。, 并 进行 了 相关 分 
析 。 其 中 ,由 联合 国 粮农 组 织 FAO(1998 年) 推荐 的 
Penman- Monteith 模型 趾 , 既 考虑 了 作物 的 生理 特 
征 , 又 考虑 了 空气 动力 学 参数 的 变化 ,具有 较 充 分 
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的 理论 依据 和 较 高 的 计算 精度 ,得 到 较为 广泛 的 应 
用 一 ”。 目 前 针对 西藏 地 区 ET 的 研究 多 停留 在 局 
部 区 域 .少量 站 点 、 较 早 时 段 ,内 容 主 要 涉及 年 E 
节 ETo 时 空 变化 趋势 及 成 因 分 析 等 方面 。 例 如 杜 军 
等 利用 西藏 “一 江 两 河 "主要 农 区 4 个 气象 站 点 资 
料 , 分 析 了 农田 ET 空间 分 布 .年 代 际 变化 特征 以 及 
讨论 了 影响 ET 变化 的 气象 因子 ; 杨 永 红 等 “重点 
分 析 了 7 个 站 点 逐日 . 逐 月 ET 年 内 变化 规律 ; 毛 飞 
等 的 分 析 了 那曲 地 区 气温 .降水 和 潜在 蒸 散 的 地 理 
分 布 以 及 年 内 ,年 际 变化 规律 ; 何 慧 根 等 “讨论 了 
动力 .热力 和 水 分 因子 对 洪 在 蒸 散 的 影响 ; 张 娜 等 ” 
采用 曼 - 肯 德尔 法 (Mann-Kendall,M-K) 对 西藏 年 际 
ET。 进 行 趋势 检验 及 空间 插值 分 析 ; 唐 小 萍 等 ”基于 
全 球 气候 变 暧 背景 ,分 析 了 西藏 雅 江 中 游 地 区 ET) 
时 空 分 布 特征 以 及 贡献 率 。 

综 上 ,已 有 的 研究 对 西藏 尤其 是 西藏 不 同 气 候 
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区 ET。 及 未 来 ET 变化 的 研究 较 少 ,本 文 根 据 温 克 刚 
等 ”的 研究 ,将 西藏 划分 为 那曲 中 西部 等 五 大 气候 
区 ,通过 分 析 各 气候 区 ET 的 变化 特征 ,探讨 其 突变 
性 周期 性 及 灰色 模型 对 未 来 ET 的 适用 性 , 旨 在 为 
西藏 节 水 灌溉 合理 调配 水 资源 及 干旱 监测 提供 科 
学 依据 。 


1 材料 与 方法 


1.1 研究 区 概况 与 数据 来 源 

西藏 位 于 祖国 西南 边 卫 , 是 青藏 高 原 的 主体 ， 
地 理 坐 标 在 26° 50’ ~36° 53’ N 78° 25’ ~99° 06'E Ż 
Ho ARAR 、 积 温 小 .昼夜 温差 大 . 干 湿 季 明 显 , 年 
平均 气温 为 -2.4~12.1 "C ,年 降水 量 66.3~894.5 mm, 
年 日 照 时 数 1443.5~3574.3 ho 

本 文 所 需 气 象 数据 来 源 于 西藏 自治 区 气象 局 
提供 的 1981—2019 年 西藏 地 区 38 个 气象 站 点 的 逐 
日 气象 资料 ,包括 平均 气温 、 最 高 气温 .最低 气温 、 
平均 风速 .日照 时 数 和 平均 相对 湿度 。 其 中 ,西藏 
根据 气候 区 域 划分 为 那曲 中 西部 (那曲 ` 班 戈 、 申 
FL GME KE OR aT .普兰 . 改 则 ) 沿江 一 线 ( 拉 
BY 泽 当 .日喀则 I IR EP EARS HR, 
SATE .加 查 PAAR) .南部 边缘 地 区 ( 聂 拉 木 、 
定 日 .隆子 、 错 那 , 帕 里 )、 东 北部 (昌都 、 丁 青 、 索 县 、 
嘉 黎 .比如 洛 隆 .类 乌 齐 ) 和 东南 部 (林芝 波 密 、 察 
隅 \ 左 页 、 芒 康 、 米 林 、 八 窒 )( 图 1)。 
1.2 研究 方法 
1.2.1 Penman-Monteith 修正 公式 洪 在 蒸 散 量 
ET 由 联合 国 粮农 组 织 FAO(1998 年 ) 推 荐 的 彭 曼 公 


>z 


那曲 中 西部 


i 和 东北 部 。 
“沿江 一 线 žá a a 


ee es wns: 
O 行政 边界 0 100km 
注 :基于 审 图 号 为 藏 S(2021)007 号 的 标准 地 图 制作 ， 
底 图 无 修改 。 下 同 。 
图 1 西藏 地 区 38 个 站 点 及 5 个 分 区 的 空间 分 布 
Fig. 1 Spatial distribution of the 38 stations and 


5 zones in Tibet 


式 计算 得 到 ,计算 公式 如 下 : 
900 
r- 0.408A(R, - G) + Ya 97g Hale, -e,) m 
° A+y(1 +0.34U,) 

IEP ET, AACE ZS GE ; Ry pe aT (MJ m 
d') GA ASE PGE (MJ -m?-d') ; 0,272 m 高 处 风 
速 (m…s ie 为 饱和 水 汽 压 (kPa) ;e, 为 实际 水 汽 压 
(kPa) ;A 为 饱和 水 汽 压 曲线 斜率 (kPa C); y 为 干 


湿 表 常数 (kPa.C-')。 各 因子 具体 计算 公式 参照 广 
a, 


1.2.2 小 波 分 析 和 突变 分 析 

本 文选 用 Morlet 小 波 分 析 ETo 在 时 频 两 域 的 变 
IERIE o Mann-Kendall 的 检验 方法 是 非 参数 统计 
方法 ,该 方法 不 易 受 异常 值 干扰 ,样本 无 需 遵从 一 
定 的 分 布 ,可 以 用 于 对 ET, 的 趋势 显著 性 分 析 和 突 
AS o 
1.2.3 自然 正 交 分 解 (EOF) ”自然 正 交 函 数 (又 称 
经 验 正 交 函数 ) ,无 固定 的 函数 形式 ,图 形 是 由 场 本 
吴 来 决定 的 , 它 具 有 收敛 快 又 能 更 好 地 反映 出 场 的 
基本 结构 特征 。EOF 分 析 既 可 以 取 空 间 不 同 站 点 
进行 分 解 ,也 可 以 对 同一 站 点 的 不 同时 间 .不 同 高 
度 等 多 种 要 素 进 行 综合 分 析 。EOF 分 解 成 的 时 间 
和 空间 部 分 ,具有 正 交 性 的 特点 ,用 前 几 个 方差 贡 
献 比 较 大 的 空间 函数 ( 模 态 ) 即 可 拟 合 出 原 矩 阵 的 
主要 特征 天 2 。 
1.2.4 灰色 模型 预测 “灰色 预测 的 基本 思路 是 通过 
对 某 一 个 可 看 作 灰 色 系 统 的 研究 对 象 ,建立 灰色 模 
型 进行 预测 。 灰 色 系 统 模型 的 建立 无 需 大 量 样本 
数据 和 有 规律 的 样本 分 布 , 且 工 作 量 小 预测 精度 高 ， 
逐渐 被 用 于 近期 .短期 以 及 中 长 期 的 系统 预测 。 
本 文 基于 Matlab 实现 灰色 预测 过 程 , 包 括 对 原始 数 
据 进 行 累 加 构造 累加 矩阵 与 常数 向 量 .求解 灰 参 数 
以 及 将 参数 代入 预测 模型 进行 数据 预测 等 内 容 *。 


2 结果 与 分 析 


2.1 ET, 时 空 分 布 特征 

2.1.1 西藏 ET 的 变化 将 征 ” 近 39a 来 西藏 那曲 中 
西部 和 沿江 一 线 ET 的 多 年 均值 比较 接近 ,分别 为 
663.25 mm 和 728.91 mm, ,南部 边缘 和 东南 部 ET 较 
接近 ,分别 为 510.14 mm 和 537.42 mm, 而 东北 部 最 
少 ,为 436.28 mm。 全 区 年 均 ET, X 597.12 mm, 
1981—2001 年 呈 显 著 的 减少 趋势 ,2002 一 2019 年 呈 
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显著 的 增加 趋势 (P<0.01) ,下 降 、 上 升 率 分 别 为 
88.41 mm .(10aj ', 109.70 mm (10a)'( 图 2a);2 个 不 
同时 间 段 的 高 值 区 均 分 布 在 那曲 中 西部 和 沿江 一 
线 , 其 余人 研究 区 属于 低 值 区 (图 2b~e)( 其 他 气候 区 
图 省 略 )。 
2.1.2 西藏 年 ET 的 了 OF 分 析 运用 经 验 正 交 函数 
(EOF) 对 西藏 地 区 年 ET, 的 标准 值 进行 分 解 , 分 析 结 
果 表 明 , 前 5 个 特征 癌 量 特征 值 的 累积 贡献 率 达 到 
81.3%, 但 只 有 前 3 个 特征 根 的 误差 范围 不 重 苹 。 通 
过 North 显著 性 检验 , 累积 贡献 率 接近 71.3% , 因此 
这 3 个 特征 根 可 以 很 好 地 解释 西藏 近 39 a 的 ETv 主 
要 时 空 分 布 特征 。 第 一 模 态 (EOF1) 特 征 向 量 的 方 
差 贡 献 率 为 51.5% ,在 整个 研究 区 域 均 为 正 值 ,表明 
1981 一 2019 年 西藏 的 年 ETv 变 化 趋势 在 空间 上 具有 
高 度 的 一 致 性 。 年 ET 高 值 中 心 位 于 沿江 一 线 , 低 
值 中 心 位 于 南部 地 区 ,整体 空间 分 布 呈 现 出 由 中 部 
中 心 向 西南 和 东南 逐步 递减 的 分 布 特征 (图 3a)。 
在 经 验 正 交 函数 (EOF) 分 解 中 ,时 间 系 数 为 正 
时 ,表现 为 西藏 地 区 年 ET, 增 加。 第 一 模 态 (EOF1 ) 
时 间 系 数 (图 3d) 可 以 看 出 ,1981 一 1989 年 .2007 年 、 
2009 一 2010 年 .2012 年 .2014 一 2019 年 的 时 间 系 数 
大 于 零 ,表明 这 5 个 时 段 西 藏 地 区 年 ET, 呈 上 升 趋 
势 ,这 与 西藏 年 蒸 散 的 变化 趋势 相 吻 合 , 反 之 ,其 他 


1981 一 2001 年 
线性 (1981 一 2001 年 ) 


800 F (a) 


= 2002—2019 
sis 线性 (2002 一 2019 年 )” a" 


a 
© 
© 

> 


a 
500 上 y=-8.841x+18191, 
R°=0,806 


400 1 1 1 1 


时 间 段 年 ET, 旱 下降 趋势。 

第 二 模 态 (EOF2) 的 贡献 率 为 12.19 ,西藏 年 
ET 在 空间 上 表现 为 西藏 中 西部 正 值 中 心 与 西藏 东 
南部 负 值 中 心 呈 相反 的 空间 变化 特征 ; 正 值 主要 分 
布 在 阿里 地 区 日喀则 中 西部 和 那曲 市 大 部 ,其余 
为 负 值 区 (图 3b)。 再 结合 EOF2 时 间 系 数 可 以 看 出 
(Al 3e) ,1981 一 1990 年 主要 表现 为 西藏 中 西部 的 
ET 减少 ,而 东南 部 增加 ,1990 一 2005 年 东西 部 ET。 
减 增 交替 变化 ,直到 21 世纪 10 年 代 中 期 开始 出 现 
相反 的 变化 趋势 , 即 西藏 中 西部 ET 大 幅度 增加 ,而 
西藏 东南 部 呈 递 减 趋势 。 

第 三 模 态 (EOF3 ) 的 贡献 率 为 6.0% ,该 模 态 下 
西藏 ET 呈现 中 南 .东南 正 值 中 心 与 西北 部 负 值 中 
心 的 分 布 特征 ; 正 值 主要 集中 在 昌都 和 林芝 市 大 
部 日喀则 市 中 部 及 山南 市 东部 ,其 余 区 域 为 负 值 
区 (图 3c)。 从 EOF3 时间 系 数 发 现 ( 图 3f) ,1981 一 
2019 年 EOF3 时 间 系 数 在 正 负 间 振 荡 , 表 明 西 藏 地 
区 ET, 年 际 间 的 变化 趋势 显著 。 

2.2 ET, 突 变 分 析 

2.2.1 BAFFET. ZLA AH ”西藏 地 区 38 个 站 
点 1981 一 2019 年 年 ET 的 Kendall 趋势 系数 空间 分 
布 如 图 4a, 从 图 中 可 以 看 出 ,Kendall 趋势 系数 为 负 
值 的 气象 站 点 有 15 个 ,其 中 贡 嘎 隆子 、 八 答 站 ET。 
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Fig. 2 Interannual and spatial variations of ET» in the whole region 
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(a) EOF1 
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图 3 西藏 年 ET 的 经 验 正 交 隐 数 (EOF ) 模 态 分 析 
Fig.3 EOF modal analysis of ET, in Tibet 


下 降 趋 势 非常 显著 ,其 Kendall 趋势 系数 分 别 为 
-4.63 ,-2.67 和 -4.53。Kendall 趋势 系数 为 正 值 的 
台 站 共 23 个 , 占 总 站 点 数量 的 61% , 主要 集中 在 那 
曲 中 西部 和 南部 边缘 地 区 , 且 西 藏 最 长 连续 无 降水 
日 数 的 高 值 区 主要 位 于 干旱 与 半 和 干旱 气候 区 品 ,可 
见 该 区 域 气候 出 现 了 暧 干 化 现象 ,这 与 王 娜 ”提出 
的 “西藏 部 分 地 区 呈现 暖 干 化 趋势 ( 藏 北 . 藏 西 北 干 
FA EPABX ) ,气温 升 高 带 来 蒸发 量 和 于 旱灾 害 大 
机 率 发 生 ” 观 点 相 吻 合 ;此 外 , 狮 录 河 、 聂 拉 木 i 
萨 、 帕 里 、 索 县 站 和 米 林 站 的 上 升 趋势 均 达到 显著 
性 水 平 (Poos=1.96)。 

2.2.2 西藏 年 ET 突变 分 析 利用 Mann-Kendall 方 
法 对 研究 区 38 个 气象 站 点 进行 突变 分 析 ,并 结合 滑 
动 1 检 验 和 Yamamoto 法 进行 校 验 (图 4b)。 本 文 针 
对 Yamamoto 法 ,基于 时 间 序 列 (1981 一 2019 年 ), 人 


为 设置 某 一 时 刻 为 基准 点 ,基准 点 前 后 长 度 分 别 取 
m=m=3 m=12=5 ,结果 发 现 仅 那曲 站 和 当 雄 站 发 生 2 
次 突变 ,20 个 站 发 生 1 次 突变 ,其 余 站 未 发 生 突变 。 
38 个 气象 站 点 的 突变 年 份 中 ( 表 略 ) ,在 20 世 纪 80 
年 代 发 生 突变 的 台 站 最 多 (13 个 ), 其 次 是 21 世纪 
10 年 代 (8 个 )。 发 生 突变 的 站 点 主要 分 布 在 南部 边 
缘 地 区 .沿江 一 线 和 东北 部 ,其 中 南部 边缘 地 区 所 
属 站 点 全 部 发 生 突变 ,生出 现年 ET, 上 升 的 突变 ( 除 
隆子 站 ) ,表明 该 区 域 年 ET, 变 化 具有 强 敏 感性 ;未 
发 生 突变 的 站 点 主要 集中 在 那曲 中 西部 。 结 合 
4a 可知, 发 生 突 变 的 Kendall 系数 正 负 值 数量 均等 ， 
表明 这 些 站 点 的 ERT, 经历 了 上 升 或 下 降 变化 。 
2.3 ET, 周 期 变化 特征 

为 进一步 了 解 西藏 地 区 年 ET 的 周期 性 ,本 文 
利用 Morlet 小 波 变 换 对 1981 一 2019 年 全 区 及 5 大 气 
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图 4 西藏 地 区 ET 的 Kendall 趋势 系数 和 气象 站 点 突变 分 布 


Fig.4 Kendall trend coefficient of ET, and mutation distribution of weather stations in Tibet 


候 区 年 ETo 时 间 序 列 进行 分 析 , 全 区 ,沿江 一 线 年 均 
ET 在 30~35 a 时 间 尺 度 振荡 较为 明显 ,表现 为 “ 正 - 
负 ” 交 蔡 ,表明 近 39 a 来 全 区 .沿江 一 线 ETv 经 历 了 
“ 偏 大 - 偏 小 "交替 过 程 ;而 在 11 a 时 间 尺 度 也 存在 
较为 密集 的 正 负 交替 ,表明 该 区 域 ET 在 短 时 间 尺 
度 上 存在 更 多 “ 偏 大 期 "和 “ 偏 小 期 "(图 5a、5c)。 通 
过 小 波 方差 检验 确定 全 区 沿江 一 线 年 均 ET 在 33a 
时 间 尺 度 上 达 最 高 小 波 方差 , 即 为 第 一 主 周期 ,此 
外 在 11a 的 时 间 尺 度 上 对 应 第 二 主 周期 (图 5g、 
5i)。 那 曲 中 西部 年 平均 ET, 在 30~35 a、20~25 a 和 
3~8 a 时 间 尺 度 呈 现 交 替 过 程 (图 $b) ,加 之 方差 验 
证 (图 5h), 第 一 主 周期 为 33 a, 第 二 .三 主 周期 分 别 
在 22a 和 3 a。 图 5d 表 明 南 部 边缘 地 区 在 30~35 a 
存在 正 负 交替 的 周期 振荡 ,而 在 8~12 a 存在 密集 正 
负 交 蔡 , 进 一 步 通过 图 5j 可 见 , 小 波 方 差 图 有 3 个 峰 
值 ,其 3 个 主 周期 分 别 是 33 a、8a 和 12a。 图 5e 中 ， 
东北 部 大 时 间 尺 度 在 30~35 a, 目前 处 于 偏 大 时 期 ， 
小 时 间 尺 度 3~8 a 较为 密集 ,结合 图 5k 来 看 ,33 a 和 
8 a 分 别 是 第 一 、 二 主 周期 。 图 5f 显 示 东 南部 大 时 间 
尺度 在 30~35 a, 小 时 间 尺 度 3~5 a、8~12 a 时 间 尺 度 
交替 相对 密集 ,通过 图 51 检 验 , 小 波 方差 图 有 4 个 峰 
值 , 共 4 个 主 周期 分 别 是 33 a、22 a、12a 和 5 as 
2.4 灰色 预测 模型 

灰色 预测 模型 CM(1, 1) 建 立 的 条 件 , 当 且 仅 当 灰 
参数 ae (-2,2) , 且 原 始 序列 的 级 比 oX k) < (0.1353, 
7.389) 全 部 满足 小 。 本 文 拟 通过 对 1981 一 2019 年 全 
区 38 个 站 点 年 均 ETo 日 期 序列 建立 GM(1, 1) 模 型 、 
级 比 和 残 差 检验 ,得 到 各 站 点 模型 中 的 参数 等 结 
果 ( 表 1)。 由 表 可 见 ,该 模型 灰 参 数 w>-0.3 且 cue 


(-2, 2) ,加 上 级 比 是 平滑 的 ,符合 灰色 系统 理论 范 
BE; 预测 精度 检验 在 80% 以 上 的 站 点 有 33 个 (25 个 
三 级 .8 个 二 级 ) , 占 研究 站 点 总 数 的 87% , 故 本 文 构 
建 的 GM(1, 1) 模 型 适用 于 西藏 年 均 ET, 日 期 序列 的 
中 长 期 预测 。 

剔除 不 合格 建 模 的 5 个 站 点 ,分 别 构建 全 区 及 5 
大 气候 区 的 灰色 模型 ( 表 略 ) ,并 根据 各 预测 模型 分 
别 得 出 西藏 2020 一 2039 年 和 2020 一 2050 年 年 均 
ET 的 预测 值 。 分 析 发 现 , 东 南部 年 均 ETo 有 明显 下 
降 趋 势 ,其 他 区 域 均 呈现 上 升 态 势 ( 南 部 边缘 地 区 
上 升 最 快 ), 且 两 组 预测 值 的 气候 倾向 率 基本 无 变 
化 。 以 上 结论 与 赵 俊 芳 等 ”在 IPCC 温 室 气 体 未 来 
排放 情景 A2( 国 内 或 区 域 资源 情景 ,其 输出 结果 作 
为 预测 未 来 2011 一 2050 年 西藏 地 区 干 湿 状 况 的 输 
入 数据 ) 下 得 出 “ 除 零 星 地 区 外 ,西藏 地 区 参考 作物 
蒸 散 量 均 呈 增 加 趋势 ,是 增长 幅度 较 大 的 地 区 位 于 
日 喀 则 市 南部 的 边缘 地 带 ” 相 吻合 。ET 的 变化 必 
然 是 受到 西藏 气候 的 影响 ,最终 其 变化 趋势 将 直接 
影响 到 当地 的 水 资源 及 生态 状况 。 


3 讨论 


本 文 基于 EOF 分 析 、 小 波 分 析 和 灰色 预测 模 
型 ,深入 研究 西藏 地 区 及 5 大 气候 区 ET, 的 时 空 演变 
规律 及 未 来 变化 趋势 。 结 果 表 明 ;西藏 地 区 近 39 a 
ET, 整 体 上 呈 上 升 趋势 ,1981 一 1988 年 处 于 ET 大 值 
区 ,1989 年 后 处 于 低 值 区 ,2005 年 起 又 持续 回升 ,这 
与 张 璐 等 ”对 锡 林 河 流域 、 称 文彬 等 "对 若尔盖 湿 
地 、 张 娜 等 ”对 西藏 ET 得 出 的 结论 较 一 致 ,与 拉巴 
等 对 藏 北 那 曲 、 钟 巧 等 "对 博 斯 腾 湖 流域 平原 区 
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图 5 西藏 ET, 在 全 区 及 5 大 气候 区 的 Morlet 小 波 变换 实 部 时 频 分 布 和 小 波 方 差 检 验 


Fig.5 Time-frequency distribution of real part of Morlet wavelet transform and wavelet variance test of ET) in the whole region 


and the five major climate regions in Tibet 


ET 得 出 的 结论 不 同 , 可 见地 域 差异 会 对 变化 趋势 
产生 一 定 的 影响 。 西 藏 中 西部 为 蒸 散 高 值 中 心 ,在 
干旱 、 半 干旱 区 ,植被 稀疏 可 能 易 导 致 强 燕 发 ;西藏 
东南 部 的 ETo 明 显 小 于 中 西部 地 区 ,该 区 域 多 为 森 
林 和 草原 区 ,具有 相对 较 好 的 植被 分 布 以 及 湿润 和 
半 湿 润 气候 ,表现 出 ETv 低 值 区 ,这 与 杨 永 红 等 "对 


西藏 高 原 区 BA SE EO h E PRL KC BLD K 
域 较 吻 合 

那曲 中 西部 在 2009 年 前 呈 缓 慢 下 降 趋 势 (图 
略 ) ,这 与 草 嘎 等 的 关 于 藏 西 北 地 区 潜在 蒸 散 呈 现 
显著 减少 趋势 的 结论 基本 吻合 ;沿江 一 线 上 升 率 极 
小 ,1981 一 2005 年 表现 为 不 同 程度 的 减 小 趋势 (图 
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表 1 各 站 点 GM(1, 1) 模 型 的 参数 和 检验 结果 
Tab.1 Parameters and test results of GM(1, 1) model at each site 


气候 区 。 站 点 模型 中 的 参数 wri arto” | 气候 区 。 站 点 模型 中 的 参数 精度/%。 ito” 
a u a u 

AS “那曲 -0.0036 424.38 86.55 0.61~1.66 | 南部 “有 聂 拉 木 -0.0076 421.78 91.60 0.78~1.21 

HH ges -0.0016 657.56 ”89.71 ”0.62~1.51 | WH 定 日 -0.0025 670.98 8846 0.80~1.31 

中 地 区 

西 扎 0.0012 710.82 87.12 0.65~1.32 隆子 0.0037 752.32 93.24 0.80~1.19 

部 ME -0.0043 443.73 81.01 0.58~1.63 错 那 -0.0064 258.23 89.19 0.76~1.31 
安 多 0.0008 508.92 87.29 0.65~1.44 帕 里 -0.0071 310.47 88.95 0.63~1.27 
狮 果 河 -0.0061 740.05 91.36 0.77~1.24 || Æ 林芝 -0.0006 483.09 87.46 0.65~1.82 
普兰 -0.0011 784.19 94.75 0.80~1.18 南 波 密 -0.0034 353.26 87.46 0.77~1.35 
改 则 0.0003 862.05 92.26 0.78~1.37 pay 0.0002 575.12 81.51 0.74~1.52 

沿 “拉萨 -0.0079 693.78 90.42 0.46~1.62 左 贡 -0.0031 403.29 86.99 0.74~1.78 

w aM 0.0047 1001.80 87.58 0.57~1.78 芒 康 0.0007 467.01 85.04 0.66~1.27 

5 日 喀 则 -0.0105 526.21 85.61 0.70~1.40 米 林 -0.0034 327.10 90.58 0.79~1.31 
AA 0.0009 688.26 86.10 0.68~1.28 八 宿 0.0073 1227.70 89.23 0.53~1.42 
浪 卡 子 -0.0009 592.48 90.30 0.70~1.34 | 东 昌都 -0.0148 349.23 74.83 0.62~1.42 
Piik -0.0003 884.31 79.57 0.46~1.62 | 北 TH -0.0055 321.21 79.52 0.68~1.84 
尼 木 -0.0060 514.40 63.90 0.60~1.87 i RH -0.0083 369.50 80.44 0.64~1.68 
Tel 0.0113 924.95 88.36 0.56~1.35 嘉 黎 -0.0038 298.42 87.83 0.73~1.39 
墨 竹 工 卡 -0.0012 792.12 88.23 0.71~1.32 比如 0.0059 427.13 83.81 0.66~1.21 
加 查 -0.0012 634.69 89.50 0.73~1.43 洛 隆 0.0022 732.76 88.66 0.66~1.36 
南 木 林 0.0063 678.85 59.13 0.44~1.44 类 乌 齐 0.0019 375.38 85.66 0.62~1.57 


注 :参数 a 为 GM(1,1) 预 测 模型 待定 系数 ;a 为 模型 的 发 展 系数 ;u 为 内 生 控 制 系数 。 


略 ), 这 与 杜 军 等 31971 一 2005 年 西藏 “一 江 两 河 ” 
结论 较为 一 致 。 文 中 EOF1 F , “1981—1989 年 ET, 
旦 增加 态势 .1990 一 2006 年 呈 减 少 趋势 "正好 印证 
了 和 毛 飞 等 I 对 那曲 市 1961—2000 年 潜在 蒸 散 
“1983—1989 年 为 偏 高 期 .1989 年 以 后 为 偏 低 期 ”的 
研究 结论 ; 相 较 唐 小 萍 等 2 研究 ,拉萨 ETo 呈 明显 夫 
加 趋势 江 孜 和 泽 当 ET 分别 表现 不 显著 和 显著 的 
减少 趋势 与 本 文 结论 吻合 ,但 日 喀 则 结论 相反 ,可 
能 由 于 研究 年 份 存在 较 大 差别 (1961 一 2009 年) 的 
缘故 。Kendall 趋势 系数 为 正 值 的 台 站 偏 多 ,与 杨 永 
红 等 得 出 的 结论 相反 , 究 其 原因 本 文 涉及 38 个 站 
点 均值 ,而 杨 永 红 研究 选择 的 是 5 个 不 同年 代 际 的 7 
个 典型 站 点 , 故 各 有 缘由 解释 ;38 个 站 中 发 生 突变 
的 时 间 集 中 在 20 世 纪 80 年 代 , 与 张 娜 等 "得 出 的 
西藏 各 气候 区 Pettitt 法 突变 结果 基本 吻合 。 

本 文 着 重 对 西藏 5 大 气候 区 ET 的 变化 特征 .其 
突变 性 .周期 性 及 未 来 ETu 的 走势 做 了 详尽 分 析 , 但 
由 于 当前 西藏 保持 多 年 连续 气象 资料 的 国家 标准 
气象 监测 站 密度 稀 玖 、 空 间 分 布 不 均 , 新 增 站 点 时 
间 序 列 较 短 等 因素 ,进行 ET 突变 和 周期 性 分 析 时 


TEK 


ZTE , ,后 续 研 究 尤 其 是 针对 不 同 区 域 时 应 合 
理 增 加 前 期 时 间 序 列 。 另 外 ,虽然 本 文 构建 的 GM 
(1, 1) 预 测 模型 大 多 通过 残 差 检验 ,但 均 未 达到 一 
级 (好 ) 等 级 ,上 且 出 现 5 个 不 合格 的 样 点 ,后 期 需 进行 
残 差 修正 ,以 提高 模型 的 适用 性 。 


4 结论 


(1) 本 文选 取 38 个 站 点 相同 时 间 序 列 最 新 数 
据 ,围绕 西藏 不 同 气候 区 进行 ET 的 相关 研究 ,从 时 
间 尺 度 来 看 ,西藏 地 区 ET, 多 年 平均 值 在 1981 一 
2001 年 时 间 段 为 588.05 mm, 下 降 速率 为 88.41 mm: 
(10a) '; 在 2002—2019 年 时 间 段 为 607.71 mm, EFF 
速率 为 109.70 mm:(10a)"', 晶 2 个 不 同时 间 段 的 高 燕 
散 区 均 集中 在 那曲 中 西部 和 沿江 一 线 。 从 空间 格 
局 来 看 ,西藏 地 区 年 ET 主要 呈现 由 中 部 中 心 向 西 
南 和 东南 逐步 递减 的 分 布 特征 ,高 值 中 心 集中 在 沿 
江 一 线 , 低 值 中 心 位 于 南部 地 区 ;其 次 在 空间 上 表 
现 为 西藏 中 西部 正 值 中 心 与 西藏 东南 部 负 值 中 心 
反 回 的 变化 特征 。 
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(2) 正 是 由 于 增加 了 2013 一 2019 年 时 间 序 列 ， 
Kendall 趋势 系数 发 生 了 较 大 变化 ,得 出 西藏 地 区 贡 
E 隆子 、 八 窒 站 ET 下降 趋势 非常 显著 (Kendall 趋 
势 系 数 <-1.96); 正 值 区 主要 集中 在 那曲 中 西部 和 南 
部 边缘 地 区 , 且 狮 泉 河 、 聂 拉 木 .拉萨 、 帕 里 、 索 县 站 
和 米 林 站 的 上 升 趋势 均 达到 显著 性 水 平 (Kendall 趋 
势 系 数 >1.96)。 那 曲 站 和 当 雄 站 潜在 蒸 散 量 在 1984 
年 和 2013 年 均 发 生 2 次 突变 ,其 余 台 站 发 生 1 次 或 
未 发 生 突变 ;发 生 突 变 的 站 点 主要 分 布 在 南部 边缘 
地 区 .沿江 一 线 和 东北 部 ,时 间 集 中 在 20 世纪 80 
年 代 。 

(3) 39 a 来 全 区 .沿江 一 线 潜在 蒸 散 量 的 变化 
主要 以 33 a 和 11 a 的 周期 震荡 为 主 ,那曲 中 西部 具 
有 33 a. 22a 和 3a 主 周期 ,南部 边缘 地 区 存在 33 a、 
8a 和 12 a 的 周期 ,东北 部 显著 周期 出 现在 33 a 和 
8 a, 东 南部 则 在 33 a、22 a、12a 和 5 a 产生 。 

(4) 首次 利用 灰色 系统 理论 预测 未 来 ETo, 预测 
结果 显示 ,东南 部 年 均 ETo 有 明显 下 降 趋势 ,其 他 区 
域 均 呈 现 上 升 态势 。 
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Variation characteristics of potential evapotranspiration and the forecast of 


grey model in Tibet 


SHI Jiqing', BIAN Duo', YANG Feiyun’, GAN Chenlong’, FAN Dongliang’ 
(1. Tibet Climate Center, Lhasa 850000, Tibet, China; 2. China Meteorological Administration Training Centre, Beijing 100081, 
China; 3. Maizhokunggar Meteological Bureau, Lhasa 850000, Tibet, China) 


Abstract: This paper aims to provide a scientific basis for water-saving irrigation, the rational allocation of water 
resources, and drought monitoring in Tibet, China. We use the Penman-Monteith model, trend analysis, Morlet 
wavelet analysis, Mann-Kendall tests, the empirical orthogonal function, and gray models to explore the spatiotemporal 
evolution of potential evapotranspiration (ET) in Tibet and the central and western regions of Nagqu as well as 
the trend of future ET, changes based on daily meteorological data collected from 1981 to 2019 at 38 
meteorological stations in the region. The results demonstrate that on the time scale, the average ET, in Tibet is 
597.12 mm. It also exhibited a significant decreasing trend from 1981 to 2001 and a significant increasing trend 
from 2002 to 2019 (P<0.01). In these two distinct time periods, high evapotranspiration was concentrated in the 
central and western areas of Nagqu and along the river line. In the first mode (EOF1), the annual ET) change 
trend in Tibet had a high degree of consistency in space, mainly exhibiting a gradually decreasing distribution 
from the center of the central area to the southwest and southeast. The high-value center is concentrated along the 
river line, and the low-value center was located in the southern region. Concerning space, the positive center of 
central and western Tibet and the negative center of southeastern Tibet are reversed. The stations with positive 
Kendall trend coefficients are mainly concentrated in the central and western regions of Nagqu and the southern 
marginal area, where the climate has been warming and drying. The sites exhibiting abrupt changes in annual ET, 
are mainly distributed in the southern marginal area along the river line and the northeast and were concentrated 
in time. In the 1980s, two mutations occurred in the potential evapotranspiration at Nagqu and Dangxiong 
stations, whereas one or no mutations occurred at each of the remaining stations. The annual average ET, both of 
Tibet as a whole and of the five major climate zones exhibited an increasing trend (except in the southeast), and 
the 33 year cyclic oscillation was the first main cycle. Additionally, there was a cycle of 11 a in the whole region 
and along the river line. The midwest of Nagqu had cycles of 22 a and 3 a, the southern marginal area had cycles 
of 8 a and 12 a, the northeast had a cycle of 8 a, and the southeast had main cycles of 22 a and 12 a. The 
constructed gray model (GM; 1, 1) used for prediction had gray parameters a> -0.3 and ae(—2, 2), and the level 
ratio was in a smooth state, meaning it can be used to predict future ET). The average prediction accuracy was 
87.85%, so we predicted the medium- and long-term average annual ET, date series for Tibet. The forecast results 
revealed that except for the significant decline in the forecast value of the annual ET, in the southeast, all other 
regions are exhibiting an upward trend, and the southern edge shows the steepest increase. 


Key words: Penman-Monteith model; evapotranspiration; spatiotemporal analysis; grey model; Tibet 


